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Comment on "Bethe Ansatz Results for the 
4/-Electron Spectra of a Degenerate Anderson 
Model " 

In a recent letter , the author calculates the density 
of states for 4/ electrons coupled to a conduction band in 
the framework of the Bethe ansatz (BA) solution for the 
degenerate Anderson model. It is claimed that the results 
qualitatively disagree with the results obtained for the 
same model but using a variational approach Even 
the high energy feature in the /-spectral function near 
the 4/- level energy e/, i.e. the "normal" ionization peak 
(NIP), is argued to be qualitatively different in the two 
approaches. In the following we point out that this is not 
the case. 

We concentrate on the U ^ oo, A^-fold degenerate 
Anderson model. As the Bethe ansatz for this model 
can yield exact results only in the large bandwidth limit 
B 3> |e/|, we confine our discussion to this limit. Zvya- 
gin first presents BA results which show that the NIP 
is shifted towards the chemical potential from ef in the 
Kondo limit n/ w 1, where n/ is the mean /-occupancy. 
He then states that the peak shifts in the opposite direc- 
tion in our variational calculation That this is not 
the case is obvious from our results in the Kondo limit 
presented in Appendix C of reference This is also 
shown very clearly in Fig. 5 of our handbook article . 
The NIP position is called £/ in Appendix C of reference 

. For large N and the Kondo limit this real quantity 
is determined by the equation 


where 


r{z)^N de\V{e)\'/{z-e). 

J-B 


(1) 

(2) 


ReT{e) is positive for negative e with |e| <C B (for a 
constant V{e) it is given by A^A ln(i?/|e|)). It therefore 
follows without actually solving the equation that in our 
variational calculation the NIP shifts towards the chemi- 
cal potential as in the BA solution. Zvyagin claims that 
the solution to this equation was found by us "in the 
complex form" Apparently he was confused by an 
unfortunate misprint in the second equality of our equa- 
tion which determines if. It should read je/j instead of 
if in the argument of the logarithm. From the first part 
of the equation presented above it should be obvious that 
this is in fact a misprint. 

We should mention that the results for the /-spectral 
function using our variational method essentially 
agree with the low temperature solution of the NCA 
equations and numerical renormalization group results 
As noted by Zvyagin [Q other low energy character- 
istics obtained within our scheme qualitatively coincide 
with the BA ones. 

In conclusion we have pointed out that in contrast to 
the claims in reference flli our variational results for the 


high energy feature in the /-spectral function for the 
degenerate Anderson model, which has been used suc- 
cessfully in the description of photoemission spectra of 
Cerium compounds, qualitatively agree with the BA re- 
sults. 
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